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ABSTRACT 

The aim of this paper is to vahdate a methodology for connecting the 
emission measure of individual solar coronal loops to the integrated emission 
measure of the entire solar corona, and using this connection to deduce the 
energetic properties of the corona, and then to show how this methodology can 
be applied to observations of solar- like stellar coronae. The solar validation is 
carried out by using spatially resolved X-ray observations of the Sun obtained 
from the Yohkoh satellite. This work is a further step in our effort to place 
the "solar-stellar connection" on a quantitative footing. In particular, we show 
how this analysis procedure can be used in the context of archival Einstein, 
ROSAT and EUVE data, as well as Chandra and XMM Newton data, as a 
complementary analysis tool to existing multi-thermal component models. 

Subject headings: coronae: Sun, stars. X-ray observations 
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1. Introduction 

One of the core questions of stellar coronal physics is what determines the intensity 
of both local heating of particular coronal structures and the overall heating of the corona 
taken as a whole. In order to answer this sort of question, one must look to both the Sun 
and other (solar-hke) stars: the Sun cannot tell us how the overall heating rate of a stellar 
corona depends on the presumed "control parameters" (viz., stellar rotation) since most of 
them are fixed, while in the absence of spatially resolved observations, stars cannot inform 
us directly about local heating processes. This point has motivated us in a series of recent 
papers to examine how one might go about comparing solar and stellar observations of 
what is presumably the same process - the "mechanical" heating of stellar surface layers to 
produce an X-ray emitting atmosphere (Peres, Orlando, Reale, Rosner, & Hudson, 1997, 
2000, Paper II; Orlando, Peres, & Reale 2000a; Paper I). The purpose of the present paper 
is to extend this work by showing how our detailed understanding of the energetics of 
individual solar coronal structures can be used to understand the energy budget of solar-like 
stellar coronae. 

The starting point of our analysis is the recognition that the X-ray emitting solar 
corona appears to be entirely formed by plasma magnetically confined in "loops" (Vaiana 
et al. 1973), and that while the emission from these loops can fluctuate substantially (e.g., 
Sheeley & Golub 1979; Shimizu & Tsuneta 1997), this emission is sufficiently steady on the 
relevant sound crossing and radiative cooling times that simple hydrostatic models suffice 
to describe them most of the time. Although the phenomenology of these coronal plasma 
structures seems to be well understood, their heating mechanism remains one of the riddles 
of coronal physics. Nevertheless, existing spatially resolved observations are able to tell us 
about the overall energetics of these structures, and to relate it to the ambient magnetic 
fields (Vaiana & Rosner, 1978; Title & Schrijver, 1998). A central role in our present 
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discussion will be played by the piece-wise integrated emission measure of individual loops, 

i-T+AT 

em(T) = / g(T') dT' , (1) 
Jt 

where q{T) is the differential emission measure, defined as 

,(T)..^(f)". (2) 

n is the plasma density, s is a coordinate along the line-of-sight to the coronal source, 
(Tmin, Tmax) is the temperature interval within a given loop over which we integrate the 
coronal emission, and T^ax is the loop temperature maximum. We shall assume (with little 
loss of generality) that the minimum temperature relevant to X-ray emitting loops (Tmin) 
is the same for all quasi-steady loops. The emission measure has been long known to be 
a highly useful observationally-constrained quantity for discussing the energetics of stellar 
atmospheres^ (e.g., Athay 1966, Jordan 1980). What makes the global emission measure 
especially useful is the fact that the total solar EM often does not change significantly over 
time scales of 30 minutes or more (Orlando et al. 2000b), with the obvious exception of 
flares and similar rapidly evolving, large transients. This observation agrees with the notion 
that the global emission measure for the coronal X-ray emitting matter, EM, is dominated 
by the steady X-ray emitting loops. 

^It is important to recognize that while the emission measure is useful for discussions of 
the overall energetics, its use is limited for discussions of detailed heating mechanisms: the 
latter type of discussion requires knowledge of the actual local plasma density, which is not 
available as long as the spatial substructuring of solar coronal loops remains observationally 
unconstrained. 
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Our plan is then as follows: following Paper I and Paper II, we will first discuss the 
emission measure of individual solar coronal loops in the specific context of Yohkoh X-ray 
images (§11), and show how these individual emission measures can be simply summed by 
taking advantage of elementary analytical properties of loop emission measures to obtain 
the overall coronal emission measure. Finally, demonstrate how properties of this integrated 
emission measure can be used to make deductions about the heating of the loop structures 
themselves (§111). This lattermost step then allows us to apply our analysis to the spatially 
unresolved observations of solar-like stellar coronae (§IV). 

2. Connecting the global coronal emission measure with the emission measure 



The standard Yohkoh data analysis yields pixel maps of temperature and of emission 
measure from the ratio of two images taken with different filters. From those maps one 
derives the global distribution of emission measure versus temperature, EM(T'): one divides 
the temperature range into bins of width AT, and sums the emission measure of all the 
pixels belonging to the same temperature bin (Papers I, II). The EM(T) distribution is 
then the emission measure of the plasma at temperature T summed over the temperature 
range AT, or (similar to Eq. 1) 



where Q{T) is defined as in Eq. (2); we shall assume that T varies in the range (Tmin, Tmox), 
where Tmin is the minimum temperature of material within loops contributing to X-ray 
emission, and Tmox is the maximum value of all values of the loop maximum temperature 



of single loop structures 




■T+AT 




(3) 



Tuax = max{Tinax} ■ 



(4) 
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It is worth noting that Q{T) is a presumably smooth, continuous, function while EM(T) is 
a discretized function in the form of a histogram. In our analysis we take AlogT constant, 
with 29 bins in the range 5.5 < logT(K) < 8, the nominal range of Yohkoh/SXT, i.e., 
AlogT = 2.5/29 ~ 0.09. Since AlogT is small, we approximate AT ^ TA log Tin 10 ^ rjT, 
with r] 0.2. In this section, we shall connect the emission measure distribution EM(T) to 
the distribution of emission measures for individual loops. 

2.1. The emission measure distribution for an individual loop 

First, let us turn to the emission measure of individual loops. We recall from earlier 
work (viz., Rosner, Tucker and Vaiana, 1978, in the following RTV; Maggio and Peres 
1996) that the functional dependence of the emission measure of individual loops, em, on 
temperature depends only on the maximum temperature of the given loops, apart from 
a (normalizing) multiplicative factor; furthermore, we recall that the emission measure of 
an individual steady loop, em(T), and the ascending part of the global emission measure, 
EM(T), of the entire solar corona both obey a simple power law temperature dependence, 
starting from the minimum temperature (Tmin), of the form 

T=^/2 . (5) 

In the case of the global emission measure, this was already recognized by Athay 
(1966; also Jordan 1980, Ventura et al. 1998), while Jordan (1980) - studying the very 
same emission measure - conjectured that EM(T) may be interpreted as due to the sum of 
contributions from many individual loops. More recently Laming, Drake and Widing (1995) 
have derived the whole solar disk using the EUV data collected by Malinovsky and Heroux 
(1973); the emission measure they derive is fully consistent with the above functional form. 
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The analysis to follow below shows not only that Jordan (1980) was correct, but that 
this synthesis procedure can be inverted in order to derive the properties of the underlying 
contributing loops from the integral (global) emission measure. Finally, we note that 
while the emission measure em(T) of an individual loop terminates at the maximum loop 
temperature Tmax, the global EM(T) continues for temperatures beyond that at which its 
maximum value is attained, with a functional form that also can be approximated by a 
power law of the form 

T-" , (6) 

where n is typically much larger than 3/2. 

Let us now consider the synthesis of individual loop emission measures. For convenience, 
we write the differential emission measure for a given loop as 

q{T) = bT^ ; (7) 

numerical solutions of the one-dimensional hydrostatic loop equations show that this power 
law scahng (with 7 1/2) holds well except for loops whose length far exceeds their scale 
height, or for loops with highly non-uniform local heat deposition. Since we shall not 
concern ourselves with these more complex circumstances (see below), we can use this 
power law relation to integrate as follows: 

rT+AT h . . 

em(r) = bT'^ dT' = ::^T^^^ ((1 + VT^^ - 1) brjT''^^ . (8) 

Since Tj^ax is the maximum loop temperature, we have from the above equation, 
em(ri„ax) = briT^.t^ = bT2,^r]T^^ = q{T^ax)riT^^. In our specific treatment, 
em(T) (xT^ ^3/2. ^j^^g 7 = 1/2, so that the differential emission measure q{T) 

for a given loop can be expressed in terms of the em(rinax), 
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2.2. The emission measure distribution of loops with the same maximum 

temperature 

As pointed out by many earlier authors, most of the X-ray emission from the solar 
corona comes from plasma that lies within a (coronal) pressure scale height of the surface; 
and since most coronal loops appear to have lengths that are equal to, or less than, the 
pressure scale height corresponding to their maximum temperature, this also means that 
the bulk of coronal X-ray emission derives from loops with size less than a pressure scale 
height. (This is not to say that loops that violate these conditions do not exist: such loops 
do exist, but wc argue that they do not make a significant contribution to the integrated 
coronal X-ray flux.) 

As a consequence, we can take advantage of an extremely useful fact about loops 
with roughly constant pressure: the one-dimensional equations governing the physics of 
hydrostatic loops of this type are invariant under changes of the loop field-line spatial scale 
(RTV; Maggio & Peres 1996), so that - aside from a possible multiplicative constant factor 
- the functional dependence of the loop emission measure distribution versus temperature, 
em(T'), does not depend on the loop length, but only on the loop maximum temperature 
Tmax- Thus, adding the loop emission measures for two or more loops with the same 
maximum temperature yields exactly the same functional form versus temperature, apart 
from the normalization. 

With this enormous simplification, we can group all loops of the same maximum 
temperature T^ax together, irrespective of their length (and plasma pressure) ; this group 
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may be called a Tmax-equivalence class of loops (obviously indexed by the maximum loop 
temperature for that equivalence class). Note that, since the product plasma pressure x 
loop length determines the coronal temperature (RTV), loops with very different pressure 
and length may belong to the same Tmax equivalence class. Thus, we can use the results 
of the previous subsection, Eqs. (8) and (9), to write the emission measure for the entire 
equivalence class of loops with a given maximum temperature T^ax, Eni(^max)) in the form 

Em(T, T^ax) = qu {T /T^^^f , (10) 

where Qm is the normalizing factor for the emission measure (the maximum value, occurring 
at T = Tjnax), and [3 is always approximately 3/2. We have also written Em(T, T^ax) to 
state explicitly its dependence on both T and Tmax- 



2.3. Summing the emission measure of all coronal loops 

We now model the emission measure of the entire corona as the sum of a multitude 
of static loops of different maximum temperature. If we define /(Tmax) dT^a^ as the 
probability that a loop has a maximum temperature T^ax, such that / / (Tmax) dT^^x = 1, 
and is the total number of loops0, then the emission measure vs. T of all the loops of 
equal maximum temperature Tmax, Em(T, Tmax), is given by 

Em(T, Tmax) = iV/(Tmax) em(Tmax) (T/T^^^f^, (11) 

i.e., it too is proportional to T'^^"^ up to the temperature maximum; the corresponding 
maximum emission measure is thus N f(T^s,^) em(Tmax)- 

^Equivalently one may formulate the entire treatment in terms of the amount of solar 
surface covered by the loops with temperature maximum Tmax and emission measure of the 
loop per unit area, but this would not allow us to disentangle the dependence on pressure 
(or length) from the area. 
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The overall emission measure of the corona is due to the sum of all equivalence classes 
of Tinax loops prcscut in the corona up to the largest coronal temperature, here for simplicity 
indicated as oo; this is simply 

POO 

EM(T) = / c^T^ax N f{T^^^) em(r^ax) (T/T^ax)^ • (12) 

JT 

Thus, the global emission measure at temperature T is the sum of the emission measure 
of all the loops whose maximum temperature T^ax > T . Fig. |^ shows an idealized example 
of such a sum for a distribution of loops with different Tmax^ we hypothesize a simplified 
coronal loop population made of six different equivalence classes and derive their EM(T) 
from Eq. |12[ Since the contribution of each Tmax loop equivalence class is oc T^, the total 
emission measure distribution is also oc T^, up to the lowest value of Tmax of all the loop 
equivalence classes. This will be proved more precisely in the following. The inset of Fig. |l] 
shows the /(^max) em(Tmax) of the six loop equivalence classes used for this example. 
More generally, the shape of the global emission measure distribution is determined by 
the relative contribution of each equivalence class of loops to the overall emission measure 
(Drake et al. 2000). For example, while the EM(T) distribution of the Sun has, typically, 
only one maximum, the EM(T) distribution of very active stars may show two maxima 
(e.g., Griffiths & Jordan 1998). We will comment on this point in the last section. 

2.4. The coronal EM(T) derived from observations of the solar corona 

Fig. H shows the EM(T) derived from a coronal observation made with Yohkoh/SXT. 
The ascending part of the curve, up to its peak, at the temperature Tp ^ 2 x 10^ K, is 
reasonably well approximated by a T^/'^ power law (cf. also Papers I, II, and Orlando et 
al. 2000b) up to the peak value, EMp [k, 3.5 x 10'^^ cm^^). This power law behavior 
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for the ascending part of the emission measure is present in most phases of the solar 
cycle; Paper 11 and Orlando et al. (2000b) show several examples of this behavior. If one 
examines the entire set of Yohkoh/SXT observations along the solar cycle, one finds that 
Tp varies only by a small factor, while the amplitude of the emission measure at T = Tp, 
EMp = EM{Tp), varies by several orders of magnitude; both these variations are larger at 
the time of solar maximum. 

For T > Tp, the emission measure declines, and again closely follows a single power law 
of the form T~°, with a ~ 3.5 (Fig. |^). This behavior holds for all cases of the quasi-steady 
solar corona that we have studied, although the power law index a varies as the cycle 
proceeds: a is larger (steeper power law) near solar minimum (cf. Paper II). The close 
resemblance between the synthetic emission measure distribution shown in Fig. |1] and the 
observationally-derived emission measure distribution shown in Fig. |^ suggests strongly 
that the observed distribution can indeed be regarded as a result of the superposition of 
quasi-static loops characterized by a distribution of maximum loop temperatures; in the 
following, we will then simply presume that this is indeed the case. 



2.5. Deriving the contribution of the various loop equivalence classes 

We now show how the contribution of each equivalenced class of loops, e.g., the 
unknown expression /(T^ax) em(Tmax) appearing in Eq. (12), can be recovered from 
an analysis of the observed EM(T): in Appendix A we show that this expression can be 
obtained by using Eq. to give 

M f(T ^ (T \ ^ EM(r^ax) d EM(r^ax) , 

N f{T^^^) em{T^^^) = (3 — — . (1, 

max '2-' max 
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Now, since 




for T <Tp 



EM(T) 



(14) 




) 



—a 



for T>Tp 



we obtain immediately 







for Tmax < Tp 



A^/(Tmax) em(Tmax) 



< 



(15) 




(a+l) 



for Tmax > 



The dashed hne in Fig. |^ shows the function A^/(Tmax)em(Tmax) derived from the emission 
measure distribution EM(T) shown in Fig. i.e. computed from Eq. |l^ with (3 = 1.5 and 
a = 3.5. Eq. |l^ and Fig. || show that f(T^a.x) em(Tjnax) is a sharply peaked function 
around Tp, that there are no (or, more realistically, just very few) loops with Tmax < Tp, 
and that the distribution itself decreases for Tmax > Tp as a power law whose index is 
smaller by one (i.e., is steeper) than that of EM(r). 

Fig. 1^ also shows the results of an analysis of the EM(T), using Eq. |1^ but taking 
into account the departure from power laws of the ascending and of the descending parts of 
EM(T), and considering the effect of error bars and the inherently non-linear characteristics 
of the analysis. We have made a Monte-Carlo calculation as follows: first for each data 
point of EM(T) we have generated a random value from a Gaussian distribution centered 
on the value itself with a given by the error bar; then applying Eq. |l^ to this new EM(T) 
distribution we obtain a /(Tmax) em(Tmax)- We have thus repeated the Monte-Carlo 
calculation 1000 times and determined, for each T^ax, the median value (marked with a 
diamond) and the bounds enclosing 68% and 90% of the distribution of /(Tmax) em(Tmax) 
values obtained with the simulation. These two bounds are shown as error bars of different 
lengths on each point in Fig. |^. The /(T^ax) em(Tmax) distribution obtained in this way 
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is remarkably close to the results of Eq. |T^, albeit the latter somehow underestimates the 
contribution of loops in the descending part, especially hot (Tmax ~ 10^ K) loops, and of 
loops with Tmax ~ 10^ K. The results well below 10^ K are consistent with no significant 
contribution from soft X-ray emitting loops, however with large error bars. 

For the analytical derivation to follow we will use the piecewise power law description 
made above, for ease of calculation and because it can help us to gain insight into the 
results. Also, we will discuss again the limits of our approximations for T^ax < Tp in the 
last section. 

As an aside, we note (from Eq. (11)) that each equivalence class of loops with the same 
maximum temperature Tmax will contribute to the total emission measure as follows: 



Emtot(Tmax) = / dT N /(Tmax) em(rmax) 



- max 



.|±|eM„(^)"° (16) 

where Tmin ~ 2 x 10^ K is the lowest coronal temperature we consider, and since T^in -C Tmax 
the term containing Tmin is negligible. 



3. Constraints on coronal heating 



The fact that EM(T) can be expressed as a sum over the emission measure distribution 
of many independent static (or quasi-static) loops has important implications on the heating 
of these structures. 
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3.1. The heating budget of loops with the same maximum temperature 

If Q(T) is the differential emission measure per unit temperature interval at temperature 
T of a whole class of loops of the same maximum temperature, Tmax, and P{T) is the 
radiative loss per unit emission measure, then the total radiative losses are given by 

R= dT Q(r) P(T) . (17) 

•^^min 

Since 

Q(T) = N f{T^^^) em{T^,^) x (^]^ ' , (18) 

'/-' max \ rr 



- max 



in analogy to what was discussed in §2, and P{T) ^ PqT ^ in the range 
2 X 10^ K < T < 10^ K, with Pq ~ 2.02 x 10"^^ erg s"^ cm^ K (from the 
MEKAL spectral model; Mewe, Lemen, & van den Oord 1986; Kaastra 1992; Mewe, 
Kaastra & Liedahl 1995 and references therein), 

^ ^ N /(T^ax) em(T^a.) r^-- ^^3-2 



where we have neglected T^^^^ since Tmin <^ T„ 



Po T-l , (19) 



These radiative losses are approximately one half of the total energy input to the loop, 
the other half being conducted towards the lower transition region, where it is radiated 
away (cf., Vesecky, Antiochos, & Underwood 1979). The total amount of heat delivered 
within loops with temperature maximum Tmax is thus given by 

KT^..) ^ ^J^^ITY) ^ ^^^--) em(T^ax) Po T^l (20) 
Therefore, using Eq. [T^ 

for Tmax < Tp 



2(a + /5)EMA/ Po /T^ax 



-(a+2) 

for Tmax > Tp 



(21) 
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Fig. ^ shows the heating budget h{T^s_^ for each equivalence class of loops characterized 
by the temperature maximum Tmax, as derived from the function N f {T^s_y)em{T^i^^) of 
Fig. I; recall that /5 ^ 3/2 and that, for the case shown in Fig. |, Tp ^ 1.7 x lO^K, 
EMp ^ 3.5 X 10'*^ cm~^ and a ~ 3.5. Thus, the total coronal loop heating rate up to a 
reference temperature t is given simply by 



Hit) ^2R = 2 f dT Q{T) P{T) 

•^^min 



f22) 



Our aim is now to evaluate this equation for the reference temperature t, and then to allow 
t to increase indefinitely in order to obtain the total coronal heating rate. On the basis of 
Eq. ^, Q{T) = Qm(T/Tp)^^^, Qm = ^M.m/tiTp and, for generality, we indicate with ^ the 
index of the power laws in the ascending and in the descending part of EM{T): 



p for T < Tp 
—a for T > Tp 



(23) 



From Eq. |2|, 

/■* EMm /•* / T \ ^'^ 

H{t) = 2 dT Q{T) P{T) =2 



Train 



^ ' ' ^' P{T)dT = 2^Po I T^-'dT 



T]Tp J Train \Tp 



r]Tp 



(24) 



where t = T/Tp, r^in = T^^jTp, t' = t/Tp, and 

r (t/Tp)'3-i - (T^in/Tp)'3-i 



2EMmPo , 

H{t) = X < 

r]Tp 



/3-1 

l-(Wrp)^"' , l-(t/Tp)-(-+i) 



P-1 



+ 



a + 1 



for t < Tp 



for t > Tp 



(25) 



which for t ^ Tp yields the total heating budget 



tot 



2EMmPo 
r]Tp 



1 — {Train /Tp 



1/3-1 



a + 1 



and for Tmin ^ Tp 



tot 



2EMa/Po /II 
r]Tp 1/3-1 a + 1 



(26) 



(27) 
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Fig. ^ shows the curve H{t) for the emission measure distribution EM{T) of Fig. 0. The 
above equation shows that the entire coronal heating budget is approximately proportional 
to the peak value of the emission measure distribution {EMp), and inversely proportional 
to the temperature of maximum emission measure (Tp), with a rather mild dependence on 
the indices of the ascending and descending power laws, P and —a, respectively, which just 
enter via the factor 

" ' (28) 



k a + 1 

the value of this factor is approximately 20 for the specific case shown in Fig. |^. 

For the specific case shown in Fig. H, Htot ~ 9.6 x 10^^ ergs s~^, i.e., approximately 
twice the total radiative losses we found in paper II for the same emission measure 
distribution; this is of course not a surprise since we built in this fact in our assumption 
regarding the evaluation of H immediately preceeding Eq. |2^. 



3.2. Average heating rate per unit volume within loops 

The total heat deposited in a set of loops characterized by maximum temperature T^ax 
can also be expressed as 

/i(Tmax) = ^Euj Lj ttj , (29) 

j 

where Ehj is the average heating rate per unit volume, assumed uniform inside the loop for 
simplicity, aj is the area covered by the two footpoints on the solar (or stellar) surface and 
Lj is the semilength of the loops with T^ax as temperature maximum. The index j runs 
over all loops with maximum temperature Tmax- 

We now consider the two scaling laws for the plasma confined inside a hydrostatic loop 
(RTV): 

T^ax = 1.4 X lO^ipiy/' (30) 
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Eh = 10^ p^/^ L"^/^ 



P 



(31) 



where p is the coronal pressure (by assumption, essentially constant within loops shorter 
than the pressure scale height), and Eh is the average coronal heating per unit volume 
inside such loops. These equations allow us to solve for the pressure in terms of the loop 
length and maximum temperature. 



P = T 



T 

J. rr 



L Vl.4 X 10= 

and, upon substituting Eq. |3^ into Eq. |31|, one obtains 

for Tmax < Tp 



(32) 



Eh = { 



(33) 



10^ 



T 



7/2 



L 



-2 



for Tm»^- > Tp . 



.1.4 X 103. 

where we have taken into account the fact that, in practice, very few loops exist with 
^max < Tp. We can now solve for L for T^ax > Tp from the above equation, and substitute 
into Eq. |29| to obtain, for Tmax > Tp 



h{T^..) = 9.9 X 10-^ Tli^Y.E'^1 ^ lO-^^T^/atE^l/; 



nl/2 



"3 ■ 



(34) 



3 3 

Since this expression does not depend on L, it therefore applies to all loops with the same 
maximum temperature Tmax- Equating the right-hand sides of this equation and Eq. pTI, for 



^max > Tp 



therefore 



2(a + /3)EMAf Po [T, 



r/(/?-l) T2 V T, 



-(a+2) 



Hj (^3 



(35) 



Hj tij 



2 X 10^ 



EMpPoa + /3 „ 



for Trnax < Tf 



T-("+'-'') for T^a. > Tf 



(36) 
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For the EM(Tinax) shown in Fig. a ~ 3.5; using the typical value /3 = 1.5, we find that 

for Tmax < Tp 



EE 



1/2 

Hj % 



(37) 



2 X 10' 



EMpPo^ 3.5 
-LP 

V 



T 



-7.25 



for Tma^ > Tp . 



Fig. 1^ shows, as an example, the expression J2j aj corresponding to the emission 
measure distributions EM(Tmax) of Fig. |. 

Now, in order to derive the volumetric heating rate Eh one needs to know the 
coefficients {aj} which, however, cannot be derived from the emission measure itself, nor 
can they be derived from stellar observations. 

Finally, we note that the power law dependence on Tmax in Eq. |3^ is rather steep. 
However, common experience with coronal loops shows that hotter loops (higher Tmax) are 
more intensely heated (higher E^), so that Eh oc T^^x with C > 0; thus we can conjecture 



that the steep power law dependence in Eq. ^ mostly refiects the dependence of aj on 
Tmax, i-e., that a-,- (Tmax) oc T~Jx^^ or even steeper. Equivalently, one may substitute Eh 
from Eq. 153 into to obtain 







10^ 



for Tmax < Tp 



T 



1.4 X 103 



7/2 



Y: j^L.a, = 9.7 X lO-^T^/^ ^ ^max > Tp . 



(38) 



This equation shows explicitly the link between the heating budget and the geometrical 
factors. It may at first sight appear that there is no dependence on parameters such as the 
plasma pressure; these are, however, involved through the scaling laws and their dependence 
on L. 



Appendix B presents an evaluation of the heating of loops with T„ 



2 X lO^K 
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in order to provide the reader with a feehng of the orders of magnitude involved in this 
discussion. 

4. Discussion and conclusions 

We have studied the piecewise-integrated emission measure versus temperature, 
EM{T), of the solar corona taken as a whole, i.e., as if it were an unresolved stellar corona. 
Using the facts that the corona is almost entirely composed of steady loops whose height 
is less than their pressure scale height and that the emission measure distribution versus 
temperature for each static loop is oc T'^ (with {3 a; 1.5), we have found that we can 
derive the distribution of the emitting loops versus temperature from the emission measure 
distribution EM{T), as well as the heating budget for loops with different maximum 
temperature. 

Rather than focusing on the energetics of individual loops, this paper studies the solar 
corona taken as a whole, in the same way that stellar observations treat stellar coronae. 
In our model of the solar corona, loops with the same maximum temperature Tmax are 
grouped together as their emission measure distribution as a function of temperature 
would be indistinguishable from that of a single loop with the same XJ-q^x and emission 
measure equal to the sum of their emission measures. We then sum all the contributions 
of such "equivalence classes" of loops over Tmax to obtain the total coronal emission 
measure distribution. We show further that the observed total coronal emission measure 
distribution versus temperature, EM{T), agrees with the results of such a model calculation. 
Furthermore, we show that if one differentiates the observed EM{T) with respect to 
temperature, one obtains the emission measure for all loops with a given maximum 
temperature Tmax present at that time in the corona. This suggests that our method may 
be directly applied to the stellar case, in which case loop structures cannot be resolved 
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individually, and it may represent an alternative to the fitting with multiple thermal 
components for those cases in which the emission measure distribution can be derived. 

In a recent work Aschwanden et al. (2000) study coronal loops observed with TRACE; 
they claim that these loops appear to be isothermal, i.e. very different from the standard 
scenario given by, for instance, Yohkoh data, of loops dominated by the balance of heating, 
radiative losses and, more important, thermal conduction. 

However, the temperature and emission determination by Aschwanden et al., as 
the analysis routinely used for TRACE data, is based on narrow band XUV filter ratio 
method: each of the narrow bands accepts few spectral lines which form at rather 
different temperatures, i.e. approximately between 10^ and 10^ K. As a consequence, the 
temperature vs. filter-band-ratio is a multi- valued function, and so the relevant temperature 
determination is not-unique. The solution chosen, for the data analysis routinely used, 
has been to force the 195/171 AA ratio to be (approximately) in the 0.9 — 1.8 x 10^ K 
temperature range based on the claim that this is the range of maximum line visibility at 
constant emission measure. As we discuss below, ignoring the role of emission measure in 
determining plasma visibility and the fact that loops can have widely different emission 
measure values, can lead to severe errors. 

Indeed, on the basis of coronal loop physics matured over decades of X-ray observations 
and modeling, it is known that hotter loops have a corresponding much higher pressure and 
higher emission measure compensating for the presumed lower visibility, as can be easily 
seen with hydrostatic loop models. Since the entire loop density distribution scales with 
the loop maximum temperature as T^^^, their emission measure measure should scale as 
'^max- Taking this feature into account along with the spectral response of TRACE, loops 
between 10^ K and 10^ K are all almost equally visible, much more visible than presumed 
in the TRACE data analysis system. The plasma detected in all this wide range is however 
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invariably assumed to be in the 0.9 — 1.8 x 10^ K range: any thermal structure along the 
loop is entirely wiped out and all the loops are claimed to be at the same temperature. 

Testa et al. (2001a, 2001b) have analyzed some coronal loops observed with TRACE, 
relaxing the above mentioned limitation of temperature and trying a novel approach to the 
TRACE data analysis. They have fitted both the brightness profiles in the 171 A and 195 
A bands, used for standard TRACE observations, along the loops and used the filter ratio 
values, after removing the background emission. They find that some loops appear to be at 
5 X 10^ K and not isothermal, while others are consistent with isothermal but at 10^ K, i.e. 
more ahke some loops observed with SOHO/CDS (Brekke et al., 1997) and S-055 Skylab 
(Foukal, 1976). 

Furthermore, spectroscopic work by Brosius et al. (1986), Raymond and Doyle (1981), 
Jordan (1980), Laming, Drake and Widing (1995) yield a coronal emission measure typically 
oc T^/^ and Priest et al. (2000) obtain a good fit of a coronal loop with a not-isothermal 
profile, assuming an energy balance between a heating term, heat conduction and radiative 
losses; Priest et al. claim "strong evidence against heating concentrated near the loop 
base", i.e. opposite to one of the key conclusions of Aschwanden et al. analysis. 

Incidentally, Aschwanden et al. (2000) do not show that the loops, claimed to be 
isothermal, are statistically important - i.e. very frequent - in the overall well proved 
scenario of hot coronal loops dominated by heat conduction and definitely not-isothermal. 

For all these reasons we have decided not to take into account possible "isothermal" 
coronal loops and consider just the well founded coronal loops scenario based on X-ray 
photometric observations and UV spectroscopic studies. 

The method presented here is substantially different from that of Maggio & Peres 
(1996), who carry out a "minimum x^" fit of stellar low-resolution X-ray spectra with 
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one- or two-loop coronal models: here instead we model the coronal emission measure 
distribution as due to a continuum of loops. 

Our results show that, aside from times when flares or other large transients dominate 
coronal emission, the temperature dependence of the emission measure distribution of the 
entire solar corona can be modeled as a power law of the form oc T^/^, up to a maximum 
temperature which typically lies between 10^ K and 2 x 10^ K. For T > Tp, we show that 
the emission measure distribution decreases as a steep power law in T. This decrease of 
EM(T) for T > Tp can be used to derive the distributions both of loops and of their heating 
as a function of Tmax, the loops' maximum temperature. Finally, we show (on the basis 
of Eq. |T3|) that there are relatively few loops with maximum temperature well below Tp; 



indeed, we show that the distribution of emitting loops, derived on the basis of Eq. 
and the temperature dependence of the heating rate (Eq. ^T]) are rather sharp functions 
of Tmax- An important (as yet unresolved) issue to what extent this sharp temperature 
dependence is real, rather than depending on (either) a somewhat too idealized model or 
the unavoidably limited instrumental capabilities of Yohkoh/SXT (in particular, its spectral 
resolution). Observations using EUV spectroscopic data (e.g., Brosius et al. 1996) show 
rather similar emission measure distributions, suggesting that the steep rise of the loop 
distribution is very likely real; nevertheless, it seems extraordinary that there are very few 
quasi-steady loops with T^ax < 1-7 x 10^ K. In any case, it appears reasonable that the 
contribution of quasi-steady loops with maximum temperature below 1 or 2 million degrees 
to EM(T) is indeed very small relative to the hotter corona, albeit not exactly zero, at the 
maximum of the solar cycle. 

The largest contribution to the coronal emission measure, to the coronal radiative losses 
and to the coronal heating is due to loops with maximum temperature around the peak 
of EM(T). The steep decrease of the solar EM(T) beyond the temperature of maximum 
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EM(T) implies that loops are fewer and fewer for higher and higher T^ax, and that less and 
less total heat is delivered in the hotter loops. This peak evolves as a function of the solar 
cycle: From Eq. [T^, Eq. and from the evolution of EM(T) during the solar cycle (Paper 
II; Orlando et al. 2000b) we can infer the evolution of the emitting loop distribution and of 
the related heating distribution. Observations show that both distributions become more 
concentrated at lower temperatures, and that their maximum value gets smaller during the 
minimum of the solar cycle; we also know that the ascending, low temperature, portion 
of EM(T) becomes slightly steeper during the solar minimum. The steepening during 
the minimum of the cycle may be due to several causes, including the low Yohkoh/SXT 
sensitivity to cooler plasma or (more likely) to the fact that most of the corona is devoid 
of loops and thus consist mostly of non-confined plasma. In spite of this, Tp varies only 
between ^ 2 x 10^ K and 7 x 10^ K during the solar cycle. 

On the other hand EM(T), Tp and, therefore, the distribution of loop emission measure 
and heating versus Tmax all appear to be subject to negligible (if any) changes on time scales 
of hours or - sometimes - days (away from flares). The steadiness of EM(T), and therefore 
of the distributions of the loops and of the heating, over a few hours (Orlando et al. 2000b) 
seems to imply an analogous steadiness of the heating budget of the entire confined solar 
corona, with only slow and gradual changes. It is not clear why this is so. One possibility is 
that this observational characteristic is just a consequence of summing randomly distributed 
contributions over the entire corona, so that it just reflects the effects of statistical temporal 
smoothing rather than true steadiness of the entire corona. Another possibility is that the 
global externally supplied energy budget of the corona is actually steady. Either possibility 
would have major implications for the heating process of the corona. 

Finally, we note that it is straightforward to see from Eq. |T3| that any flattening, 
or even increase, at a given temperature T* of the emission measure EM(T) along its 
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descending portion would provide evidence for another preeminent set of loops, and related 
heating, with T^ax ~ T*. In particular, a doubly-peaked EM(T) (as the one found, for 
instance, by Griffith & Jordan 1998) implies a doubly peaked loop distribution. In the 
light of our results, the findings of Griffiths & Jordan (1998) are indicative of two distinct 
distributions of steady loops peaking at two rather different temperatures, one at several 
10^ K and the other at a few 10'' K. Alternatively, it may be that steady coronal emission 
is accompanied in some (very active) stars by a relatively large number of simultanteous 
flares whose light curves overlap in such a way that they mask much - if not all - of the 
variabihty (e.g., Giampapa et al. 1996; Paper II); see Reale et al. (2000) and Peres (2000). 

Our method allows us - in principle - to infer the global distribution of the heating 
versus the loop maximum temperature; this information can, in turn, provide constraints 
on the heating mechanism of the corona and on the global coronal energy budget. On 
the practical level, however, it is far from easy to derive a relationship between the loops' 
average volumetric heating rate and the loop temperature as one needs to know the fraction 
of the solar surface covered by the footpoints of the loops with given maximum temperature 
T 

max- 

As for future applications, we plan to apply extensively the diagnostics presented here 
to the EM(T) of the whole solar corona at several points during the solar cycle, using 
the large data set collected with Yohkoh/SXT. It is also possible to apply the diagnostics 
to selected portions of the corona, e.g., Yohkoh/SXT data on single active regions, or 
to SERTS observations (Brosius et al. 1996). The diagnostics can also be modified very 
easily for application to the diffential emission measure distribution, rather than to the 
piecewise-integrated emission measure treated here. For obvious reasons, our analysis is well 
suited to stellar observations since it uses a global coronal characteristics, EM(T), readily 
obtained by a spectroscopic coronal observation. We foresee an immediate application to 
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EUVE stellar observations (Drake et al. 2000); and we also plan to apply it to data obtained 
from the spectroscopic X-ray instruments onboard the Chandra and XMM-Newton X-ray 
astronomy missions. It will be particularly instructive to compare the results for stars of 
different activity levels with the analogous solar results. 

We acknowledge useful suggestions and comments from S. Serio and E. Franciosini, 
and useful discussions with C. Kankelborg. We also acknowledge useful suggestions from an 
anonymous referee. This work has been partially supported by the Agenzia Spaziale Italiana 
and by the Italian Ministero dell'Universita e della Ricerca Scientifica e Tecnologica. 



Our task is to derive N /(Tmax) em(Tmax) from the observed EM(T). On the basis of 



A. Deriving the loop population from EM(T) 



Eq. (12), 
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since the last term is 
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and therefore 



This expression connects the distribution function of loop maximum temperatures to the 
global emission measure distribution. In the specific case of a /5 = 3/2 power law, we obtain 
instead of Eq. ^l] 



A. /(T) em(T) = ! 5M(I) _ IMEl . (A2) 



B. Example: the heating budget of loops with T^ax ~ 2 x lO^K 

In this appendix, we evaluate the heating budget of loops with Tmax in the range of 
2 X lO^K, using a few simple assumptions; our aim is to provide the reader with a feeling of 
the orders of magnitude involved. 

As a first simplification, we note that since loops with T^ax ~ 2 x lO^K largely 
dominate the emission measure distribution, their total emission measure e is comparable 
to the maximum value of EM{T) and therefore we shall assume e ~ 10^^ cm~'^. The 
emission measure in an equivalence class of loops with maximum temperature Tmax 
(all with the same length L and, therefore, pressure p) can then be approximated as 
e = UpV = p'^V/ {2kBTj^axY; rip is the proton number density, V is the loops' total 
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volume, and Ub is the Boltzmann constant; according to the RTV loop scahng laws, 
Tl^^ = (1.4 X lO-'^)''^ p L, and therefore = T^J{hL'^), where 6 « 5.7 x 10"^^ Thus, upon 
substitution, we find that the emission measure is given by the expression e = T^^^/ ipL?') 
and y = e 6 L^/T^^. 

However, loops with the same temperature may have difi^erent, even very different, 
lengths. Let us consider, then, for simplicity just two values of L, 10^ cm and 10^° cm, 
the first representative of active regions loops, the other of more quiet and extended loops 
but still shorter than the pressure scale height; their total emission measure is ei + €2 = e 
and their total volume \sV — Vi-\-V2 — h {EMiLj + EM2Ll)/T^^ (where, from here on, 
all quantities with suffixes 1 and 2 pertain to, respectively, the first and the second set of 
loops). Then if we assume that Vi = V2, and since Li = O.IL2, we are led to the result 
ei = 100 62. If we insert numerical values for the various known quantities, we then obtain 
V 10^0 cm^ 

Furthermore we may assume that loops have an aspect ratio L/r ~ 10 (a typical value 
of most loops observed) where r is the cross-sectional radius of the loop. The volume of 
each loop is then v = = 2 x ttO.I^L^ x L = 0.027rL^. On the basis of these 

considerations wc would then have rii loops for the first set and n2 for the second, such that 
ni — V/ 2vi ^ 8 X 10^ loops and 712 ~ 8 loops. 

Under the same assumptions, the sum in Eq. (38) can be re-written in the form 

^ttj/Lj = EttO.OILj , (Bl) 

so that the contribution of all loops with the same length within the T^ax 
equivalence class would amount to nn x O.OIL. Thus, for the first of loops 
considered above, riiir x O.OlLi = 8 x lO^vr x 0.01 x 10^ = 2.5 x 10^^ and 
712 TT X 0.01 L2 = 87r X 0.01 X 10^° = 2.5 x 10^. Evidently under these assumptions 
the set 1 of loops, i.e., those with L — lO^cm, have a significantly larger sum, and therefore 
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dominate the heating budget of the T^ax equivalence class. Thus, from Eq. (38) we have 
h{2 X lO^K) ^ 10-6 X 11.3 X 10^1 x 2.5 x 10"er^ ^ 28.3 x 10^^ erg. 
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Fig. 1. — A simple example of summing the piecewise integrated emission measure 
distribution of various loop equivalence classes (six in this case), of different temperature 
maxima. The thin lines yield the T^/^ contribution of each class of loops having the same 
Tmax and the dotted vertical lines mark the related Tmax values. The solid histogram is the 
sum of all the equivalence classes and is oc T^/'^ up to the lowest T^ax- The inset shows the 
assumed N f{T^^^)em{T^^^). 

Fig. 2. — Piece-wise integrated emission measure distribution for the whole corona, observed 
with Yohkoh/SXT on 6 Jan 92, at 21:45 UT (top), close to the maximum of the solar cycle. 

Fig. 3. — A^/(71nax)ewi(Tmax) derived from the EM(T) of Fig. ^. Diamonds, error bars and 
solid connecting line are the results of the Monte-Carlo calculation discussed in the text: 
diamonds mark the average values, small error tics on the error bars enclose the 68% of 
the values, the large ones the 90%. The dashed line shows the result in the alternative 
calculation, approximating the ascending part up to the maximum with a power law 
and the descending part with a T^^-^ power law, and then applying Eq. |l^ to the two power 
laws. 

Fig. 4. — The function /i(Tmax), i-e. the heating budget of each equivalence class 
corresponding to maximum loop temperature Tmax- 

Fig. 5. — The function H(t), i.e., the heating budget of the whole corona, integrated up 
to temperature t. The figure shows that most of the contribution occurs for t < Tp; Tp is 
where the slope visibly changes. 

Fig. 6. — E]{J aj for loops with maximum temperature Tmax, derived from the EM(T) of 
Fig. |. 
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Simulated EM(T) Distribution 
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Observed EM(T) Distribution 
Yohlcoh/SXT, 6 Jan 92, 21:45 UT 
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